Modified Gauss-Bonnet theory as gravitational alternative for dark energy 
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We suggest the modified gravity where some arbitrary function of Gauss-Bonnet (GB) term is 
added to Einstein action as gravitational dark energy. It is shown that such theory may pass solar 
system tests. It is demonstrated that modified GB gravity may describe the most interesting features 
of late-time cosmology: the transition from deceleration to acceleration, crossing the phantom divide, 
current acceleration with effective (cosmological constant, quintessence or phantom) equation of 
state of the universe. 
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1. The explanation of the current acceleration of 
the universe (dark energy problem) remains to be a 
challenge for theoretical physics. Among the number 
of the approaches to dark energy, the very interest- 
ing one is related with the modifications of gravity 
at large distances. For instance, adding 1/R term 
0, to Einstein action leads to gravitational al- 
ternative for dark energy where late-time accelera- 
tion is caused by the universe expansion. Unfortu- 
nately, such 1/R gravity contains the instabilities 
of gravitationally bound objects. These instabilities 
may disappear with the account of higher deriva- 
tive terms leading to consistent modified gravity^. 
Another proposals for modified gravity suggest InR 
or Trl/R terms[3L account of inverse powers of 
Riemann invariants |jj or some other modifications 
0. The one-loop quantization of general f(R) in de 
Sitter space is also done 0- In addition to the sta- 
bility condition which significally restricts the possi- 
ble form of f(R) gravity, another restriction comes 
from the study of its newtonian limit [Icj . Passing 
these two solar system tests leads to necessity of fine- 
tuning of the form and coefficients in f{R) action, 
like in consistent modified gravity That is why 
it has been even suggested to consider such alter- 
native gravities in Palatini formulation (for recent 
discussion and list of references, see [llj). 

In the present paper we suggest new class of mod- 
ified gravity, where Einstein action is modified by 
the function f(G), G being Gauss-Bonnet (GB) in- 
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variant. It is known that G is topological invari- 
ant in four dimensions while it may lead to number 
of interesting cosmological effects in higher dimen- 
sional brane- world approach (for review, see |12|). It 
naturally appears in the low energy effective action 
from string/M-theory (for recent discussion of late- 
time cosmology in stringy gravity with GB term, see 
13]). As we demonstrate below, modified f{G) grav- 
ity passes solar system tests for reasonable choice of 
the function /. Moreover, it is shown that such mod- 
ified GB gravity may describe late-time (effective 
quintessence, phantom or cosmological constant) ac- 
celeration of the universe. For quite large class 
of functions / it is possible to describe the transi- 
tion from deceleration to acceleration or from non- 
phantom phase to phantom phase in the late uni- 
verse within such theory. Thus, modified GB gravity 
represents quite interesting gravitational alternative 
for dark energy with more freedom if compare with 
f(R) gravity. 

2. Let us start from the following action: 

S = J d*xV=g(^R + f(G()) ■ (1) 

Here G is the GB invariant: G — R 2 — &R llu R tlv + 
Riiv£oR^ a ■ By introducing two auxilliary fields A 
and B, one may rewrite the action (JTJ as 

S = J d 4 x^y^R + B(G-A) + f(A)^j . 

(2) 

Varying over B, it follows A — G. Using it in (J2J, the 
action (JIJ is recovered. On the other hand, by the 



2 



variation over A in one gets B = /'(A). Hence, 

5 = | d 4 xV^ + f{A)G - Af'(A) + f(A? 

The scalar is not dynamical, it has no kinetic term 
and is introduced for simplicity. Varying over A, the 
relation A = G is obtained again. 

The spatially-flat FRW universe metric is chosen 

as 

3 

ds 2 = -dt 2 + a(t) 2 ^2(dx l f . (4) 

i=l 

The first FRW equation has the following form: 
= -4# 2 + Af'(A) - f(A) - 24Af"(A)H 3 . (5) 

K 

Here the Hubble rate H is defined by H = a /a. For 
(0J, GB invariant G ( A) has the following form: 

G = A = 24 (^HH 2 + H 4 ^ . (6) 

In general, Eq.JSJ has deSitter universe solution 
where H and therefore A — G are constants. If 
H = Hq with constant Hq, Eq.J^J) looks as: 

= -^H 2 + 24H V (24ff 4 ) - / (24tf 4 ) . (7) 

For large number of choices of the function /, Eq.Q 
has a non-trivial (Hq ^ 0) real solution for Hq (de- 
Sitter universe). Hence, such deSitter solution may 
be applied for description of the early-time inflation- 
ary as well as late-time accelerating universe. 

Let us check now how modified GB gravity passes 
the solar system tests. The GB correction to the 
Newton law may be found from the coupling matter 
to the action Varying over g pu , we obtain 

= (- R ^ + \^ R ) + T " V + l^fiG) 
-2f'(G)RR» v + Af'(G)R' J ' p R 1/p 
~2f'{G)R^ aT R» paT - Af'(G)R^ u R pa 

+2 (VV/' (G)) R ~ 2sT (V 2 /'(G)) R 
-4 (V P V^/'(G)) R vp ~ 4 (VpV/'(G)) R w 
+4 (V 2 /'(G)) R» v + IgT (V P V CT /'(G)) 
-4(V p V CT /'(G))i?^ CT . (8) 

Here T^ v is matter EMT. In the expression ©, the 
third derivative /'"(G) is included as V 2 /'(G) = 



/"'(G)V 2 G + /"(G)V^GV AI G, for example. In the 
equation JSJ corresponding to the first FRW equa- 
tion, however, the terms including /'" do not ap- 
pear. When T pv = 0, the (t, i)-component of Eq.® 
reproduces Eq. (jSJ) by identifying G with A. The per- 
turbation around the deSitter background which is 
a solution of (7J may be easily constructed. We now 
write the deSitter space metric as <7(ow 5 which gives 
the following Riemann tensor: 

R(0)^vp<T = H (g(0) W 5(0><x - 9{Q)iJ.a-9{0)vp) • (9) 

The flat background corresponds to the limit of 
Hq — ► 0. Represent g^ — g(o)p V + h pu . For 
simplicity, the following gauge condition is chosen: 

sfo'V = vf 0) v = °- Then E q-@ g ives 

= ^2 ( V ' V - 2H o V) + t m* • ( 10 ) 

The GB term contribution does not appear except 
the length parameter 1/Hq of the deSitter space 
which is determined with the account of GB term. 
Eg. (|1U|I proves that there is no correction to New- 
ton law in deSitter and even in the flat background 
corresponding to Hq — > whatever is the form of 
/. We should note that the expression i|l(J[l could be 
valid only in the deSitter background. In more gen- 
eral FRW universe, there could appear the correc- 
tions coming from /(G) term. We should also note 
that in deriving (|10f> . we have used a gauge condition 
9(0) = but if we include the mode correspond- 
ing to g? Khfj,v, there might appear corrections from 
/(G) term. Ea.lJTO)) only shows that for the mode 
corresponding to the usual graviton, any correction 
coming from /(G) does not appear in the deSitter 
background. 

In case of /(i?)-gravity 0, 0, for most interest- 
ing choices of the function / an instability was ob- 
served in 0. The instability is generated since the 
FRW equation, contains derivatives of the fourth or- 
der and therefore the scalar curvature propagates. 
This may cause the appearence of growing force be- 
tween the galaxies. Only special forms of f(R) may 
be free of such instability [j]. 

In /(G)-gravity case the scalar field denoted 
by A in © has no kinetic term. Then the scalar 
field or the curvature itself does not propagate and 
therefore there is no such instability as in |3j , which 
is clear from the (no) correction to the Newton law 
in JSJ. Thus, modified GB gravity may pass the solar 
system tests (at least, for some functions /(G)). 

Having in mind the fundamental property of the 
current accelerating universe, it is interesting to 



study the possible transition between deceleration 
and acceleration of the universe for the action (JTJ. 
First, Eq.@) could be rewritten as G = A = 
2AH 2 'd/a. Then at the transition point a = 0, the 
Gauss-Bonnet term vanishes: G = A = 0. Let us 
assume the transition occurs at t — to- The Hubble 
rate H may be expanded as 

H = Ho + H^t-toj+Hiit-to) 2 

+H 3 (t-t f + o((t-t ) 4 ) ■ (11) 

At the transition point t = to, one finds 

G = ^H 2 (H 2 +H 1 ) . (12) 
Since G = there, it follows 

Hi = -H 2 . (13) 

Hence, G (A) is: 

G = 96H 2 (H 2 - Hi) (t - t ) 

+A8H 2 (5# 4 - 2H Q H 2 + 3H 3 ) (t - t Q ) 2 
+o((t-t Q f) . (14) 

We now also assume that f(A) could be expanded 
as 

f(A) = f Q + f 1 A + f 2 A 2 + f 3 A 3 + 0(A 4 ) . (15) 

Then substituting (HJ, (T3J|, and JT^J into (JSJ, we 
find 
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48ff 2 



(16) 



Combining (|13[1 and (|16|l . one can show that the 
Hubble rate can be determined consistently, which 
suggests the existence of the transition between de- 
celeration and acceleration of the universe. Wc 
should note that H could be determined by a 
proper initial condition. Then the transition con- 
dition could be /2 7^ 0, that is, f(A) contains the 
quadratic term on A. 

Let us consider now the possible transition be- 
tween non-phantom phase and phantom phase of 
the universe (if current universe is phantom one). 
We now assume that the transition occurs at t = tx. 



Since H = at the transition point, it is natural to 
assume the Hubble rate behaves as 

H = Ho + H 1 (t-t 1 f + 0^-^ . (17) 
Hence, A (G) behaves as 

G = A = 2AH$ + 48i/ 2 i/i (t -h)+0 ((t - h) 
Eq. JSJ shows 



(18) 



T H 2 + 247J 4 /' (24if 4 ) - / (24i/ 4 ) 
1152J7J/" (24i/ 4 



Hx = 



(19) 

if /" fom^ ^ 0. We also note that H and Hi 

should be positive. Ea. ljT§|l suggests the existence of 
the consistent solution. Then the crossing of phan- 
tom divide could occur for large class of functions 
f(G). The condition could be /" ^ and H t > 
when H = Ho > 0. 

The transition between non-phantom phase and 
phantom phase could be regarded as a perturbation 
from the deSitter solution in (JJJ . The following per- 
turbation may be suggested 



H = H + SH , 
where H satisfies Q. Eq. gives 



(20) 



= SH + 3H SH 



-4JT, 



1 



+ 96« 2 ff 4 /" {2AHo l )) 5H - {2l) 



Here, it is supposed that /" (241$) ^ 0. Let A 
satisfies 



= A 2 + 3H \ - 4# 2 



1 



96k 2 H 4 /" (24ff 4 ; 



(22) 



Then 5H behaves as SH ~ e At . If A is real, 6H is 
monotonically increasing or decreasing, H does not 
vanish and therefore there is no transition. If 



600k* H$f" (24Hq) < 1 



(23) 



A becomes imaginary and SH oscillates. Therefore 
the transition between the phantom phase, where 
H = SH > 0, and the non-phantom phase, where 
H = SH < could be repeated in oscillation regime. 
The amplitude of the oscillation of SH decreases as 
\SH\ ~ c~ 3Hot / 2 and the universe asymptotically 
goes to deSitter space. 
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3. To show that modified GB gravity may lead 
to quite rich and realistic cosmological dynamics 
we consider some explicit examples of the function 
f(G). The following solvable model, which does not 
belong to the above class, may be discussed 



f(A) = fo \A\ 



1/2 



Here fo is a constant. Assuming 



H 



ho 
t -t 



(when ho > 0) 
(when ho < 0) 



Eq.JSJ gives 



= --h 2 a - 



fo (1 + hp) 
2 (h - 1) 



24/ij] fa - 1) 



1/2 



(24) 



(25) 



(26) 



The solution differs for ho > 1 and ho < cases. The 
simple analysis shows that when f§ > 3 /2k 4 , there 
are two negative solutions, which decribe effective 
phantom universe with w < -1. When / 2 < 3/2/c 4 , 
we have one ho > 1 solution, which describes the 
effective quinteessence — f < w < —1/3 and one 
ho < solution describing effective phantom. If 
16/o/3k 4 < 1/k 8 , there are two real solutions which 
satisfy < h < 1. 

For the second model let f(A) behaves as 



(27) 



in a proper limit with constants fo and /3. We now 
assume in a limit t — > oo, ff behaves as 



(28) 



with constants ho and a. If a > —1, A (G) behaves 
as 



A = G ~ 24ftft 4 



(29) 



Then © gives 4a/3 = 2a, i.e. /3 = 1/2 or a = 0. 
The case = 1/2 corresponds to the model l|24|l . On 
the other hand, since H is a constant when a = 0, 
this case corresponds to Q. Hence, the only new 
case is a < —1, where 



A = G - afrgf 



3 + 3a-l 



(30) 



Then in (J5J, the behavior of each term can be 
described as H 2 - 0(t 2a ), Af'(A) - f(A) - 
O (t( 3a -W), Af"(A)H 3 ~ O (t Q ( 3 ' 3 + 1 )). 

With the assumption 2a = (3a — l)/3 > 
a (3/3+1), one finds /3 = 2a/ (3a — 1) and a < 



— 1/3, which is consistent with a < 
sumption a < — 1 and /3 = 2a/ (3a 
1/2 < p < 2/3. Eq.® also gives 



0~% + V - 



l)/o \ahl\ 



-1. The as- 
- 1) restricts 



(31) 



which leads to non-trivial solution for h if fo < 0. 

Let us consider t — »0or£o — t — *0 limit. If 
a > —1, we obtain Ea. (|30|l . If one requires 2a = 
(3a — l)f3 = a (3/3+1), then j3 = —a and a = 
0, -1/3. In case = -a = 1/3, Eq.© tells 



°-4 + tA 



1/3 



(32) 



which requires fo > and gives a non-trivial solution 
for /i . 

If one restricts 2a = a (3/3 + 1) < (3a — l)/3, it 
follows /3 = — 1/3. In this case Eq.® gives 



3 a 32(3a-l)/ , 

= — ? /in |a/i, 

k z 3a 



3I-V3 





(33) 



where a can be arbitrary as long as a > — 1. 
If a < -1 (see Ea.ffltyl. Eq.JSJ) gives 







1 1 

2 + 4^ 
3 



24(4/3 -5)/3(/3-l) 



hi 



2(2/3-1) 



(34) 



The first equation shows 1/2 > j3 > 1/4. The second 
equation in (|34() gives a non-trivial solution for ho if 
/o>0. 

The above results show that the asymptotic solu- 
tion behaving as (|28|l with a ^ can be obtained 
only for the cases where /3 = 1/2 or 1/2 < /3 < 2/3 
when t — ► oo or /3 = —1/3 or 1/4 < /3 < 1/2 when 
i or to — t — > 0. In other situations, the asymptotic 
solution corresponds to the deSitter space Q. The 
following model may be taken as next example: 



f(A) = -- + (3A 2 



Then Eq.JTJ) gives 



0=--#o 

K 



\2Hl 



576/3#= 



(35) 



(36) 



It could be difficult to solve the above equation ex- 
actly, but if the curvature and therefore H is small, 
one gets 



HZ 



(37) 
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On the other hand if the curvature is large, we obtain 



H 



H? 



192/3k 2 



(38) 



For the consistency Hi 3> H s , hence a(3n A <C 1. 
Then the large curvature solution (|38|) might cor- 
respond to the inflation in the early universe and 
the small curvature one in l|37(l to the late time ac- 
celeration. Similarly, one can construct many more 
explicit examples of modified GB gravity which com- 
plies with solar system tests (at least, for some back- 
grounds) and leads to accelerating (effective cosmo- 
logical constant, phantom or quintessence) late-time 
cosmology. 



4. In summary, it is shown that modified GB grav- 
ity passes two fundamental solar system tests in the 
same sense as General Relativity. It easily describes 
the late-time acceleration of the universe whatever is 
the effective equation of state ( effective cosmological 
constant, quintessence or phantom). Moreover, we 
demonstrated that such modified GB gravity may 
describe the transition from deceleration to accel- 
eration as well as crossing of phantom divide. It 
remains to understand whether any limitations to 
the function /(G) from solar system gravitational 
physics may be found. More serious limitations to 
its form may be searched fitting the theory against 
the observational data. 
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